ss 


7.9 kS ate (6 
f[[5t C e 


— 


Alternatively, we can represent 








| Jw 
the complex number s in polar form s = |s| e“. 
| pE. 
> 941/2 . i 
where |s = (a — w^) ~ is the magnitude and | 
o=/s atan(w/o) the phase of s. : 
erren 

It is straightforward to derive | U 
Lf HM 2 l E 

- = —e9u|-|-- and L= = —18. 2 MES 

S  |s| S |s 5 











An advantage of laplace transform:- 


We can transform an ordinary differential equation (ODE) into an algebraic equation (AE). 


t-domain s-domain 


` . 


| © 


E” Partial fraction 
: (3) expansion 


Examples of Laplace transform 


= Unit step function nn 
f(t) = us(t) = : = O1 —— 
——— ff 


0 
^ ma —$ 1 —si| + 1 
mum F(s) = L 1-e "dt = a le w N A (Memorize this!) 


pa 


Solution to ODE 





| 
J 
m 


« Unit ramp function in: 
Je EU 
so=(0 t«0 + 
= n? aS aep VE aude i niae , i see A 
^ EAS FG)-Í[ te i^ e ¡te de Te c stas Al 
E: (Integration by parts) 
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f(t) Width = 0 


» Unit impulse function Height = inf 


Area = 1 
x eS Rz] | co l i p 
EN = f i j : f Be uc = g(05 ccc | y 
aL), | O, t = O 1 Pas NN y) g | e a e — N 1 


i 0 
>  F(s)-— l et e *9 —1 (Memorize this!) 


a , f(t) 
= Exponential function ) 
u gt 120 sod 
fe) = | 0 t«0 0 
(o) — [^ rat n SS 
MMC a F(s) = j e ted = 7. c i sa) 


« Sine function 





Remark: 
Instead of computing Laplace for each function, and/or transform for each function, 
memorizing complicated Laplace transform, use the Laplace transform table 





"Signa — Wavcformf() Transform F(s) 
Impulsc à(t) 1 
Step function u(t) > 
Ramp tu(t) 4 
Exponential [e7?' Juc(t) M 
Sta 
Damped ramp trer ]u{t) > Ll 
Sinc [sin Br]«(t) app 
c | R 
Cosine [cos Bt ]u (t) x5 
Damped sinc [e ^" sin Br luL) B _z_ 
(s + a)? + 8? 
Damped cosine [e7*' cos Br ]u(t) (s + a) 


(s + a)? + 8? 






i iil 
———$ 
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Ex: Newton's law 


de LZY f(t) 
MS = f) LM | 
— u(t) 


We want to know the trajectory of x(t). By Laplace transform, 
M (s? X (s) — sz(0) — z'(0)) = F(s) 


=> Xi = apa 2 E 


A O. 


(Total response) z (Forced response) + (Initia! condition response) 


p | A 


c LD =C O) + z(O)us(£) +2 (0)u5(t) 











The Laplace transforms of some common functions of time are shown below 





2 —üa 
te "ug(r) > 





ug(t) & 17s (say 
to 1/s? mca n! 
fe ut) < l 
i , ntl 
; (s 4 a) 
n! 
A mri o 
sint Haf e 5 


* a 1 
S(t) 1 Sro 


ó(t—a)e a” 





COSCf ugt => 


: s + 
e " ug(£) & — 
s+a -af . C 
e SIn@t D = ———— 


(sta) + (0 





-at * l 
te — ug(f) e» : 


qe f : E 
2 (s+a) e coset tolt) => Ez 


2 2 
(Stay + a 





BEA en ee Ree ee e 
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Laplace transform theorem 


Oe ere Gs 





- 
—P AS 





ltem " o Theorem Naine l 
R 
l ZM < Fo) = C flea Definition 
2 JM) < Li Linearity theorem 
3 Lf (+f, (0) = Fils) + £s Linearity theorem 
i Pef] < Pita Frequency shift theorem 
T[ft-T] ze" Time shift theorem 
Ç Tlf) = = -f Scaling theorem 
df | 
T. L E | = SP(5) - f(0=) Differentiation theorem 
g e = CBI - sf(0-) - PIL Differentiation theorem 
0 f E = 9 F(s) Leo, Differentiation theorem 
R slw =A Integration theorem 
11. flc) - ln sF(s) Final value theorem! 
12, f(0+) = lim P Initial value theorem? 
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MATLAB 


You can compute Laplace transform using the symbolic toolbox of MATLAB. If vou want to compute the 
Laplace transform of any function, you can use the following MATLAB program. 

fit) 5t, f2(t)= exp(atisii2tw) , fs(t)- 2:exp(-5t! 
clc 480 
Syms wa t = 
f1-4*t^5; a > 
disp('fl(t)s') 
pretty(f1) 
Gisp('F1(s)=") 
Fl-laplace(f1); 
pretty(F1) 
Sprintfil'in') 


f2= exp(-a*t)*sin(3*w*t); 3 Y 
disp('f2(t)=") F2(s)= | ------------—---— 
pretty(f2) 2 
disp('F2(s)=') 

F2=laplace (f2); 

pretty(F2) 
sorint£f('Win') 


£3-3*t*exp(-5*t); 
diso('£3(t}=") 


pretty (f3) 3 
disp('F3(s)=') Pg Dai mœ raeas 
F3=laplace (f3); 2 
pretty(F3) 

(s + 5) 


EX: Calculating the Laplace F(s) transform of a function f(t) is quite simple in Matlab. 


fGa)--1.25 3.5te ^" e 1.25e7* 


DB*t"espi-2"t)wil om BEG L WE Â Z 
C£, > 


— 


—-5/4/s5-+7/2/ (1s+2)"2+5/9/ (s+2) 
>> simplifylF) 

ans = 

(s-5) /s/ (st+2)%*2 


>> pretty (ans) 


which corresponds to F(s). 


2.6675). 
is um sís + 23* 


Alternatively. onc can write the function ft) directly as part of the laplace command: 


L—F2-laplace(-1.25«3.5*t*expí(-2*t)«1.25*exp(-2*t)) 







PP o Po ÓÀ— M — À— — 
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2- Inverse Laplace Transform 


The inverse Laplace transform can be obtained by use of the inversion integral given by the 
following equation 


oy 


c+ jo 
$7[F(5)] = f(t) = = f F(sje" ds,  fort>0 


The inversion integral is complicated and therefore, its use not recommended for finding 
inverse Laplace transforms To find the inverse Laplace transform of a complicated function, 
we can convert the function to a sum of simpler terms for which we know the Laplace 
transform of each term. The result is called & partial-fraction expansion. 


Partial-Fraction Expansion 


N(s 
f F (s ) = 22 Where N(s) and D(s) are polynomials in s. 
Note 


I- If the order of N(s) is less than the order of D(s), then a partial-fraction expansion can be 
made. 


I- If the order of N(s) is greater than or equal to the order of D(s), then N(s) must be divided 
by D(s) successively until the result has a remainder whose numerator is of order less than its 
denominator. 


Case 1: Partial- fraction expansion when F(s) involves distinct poles only. 
Ex: Find the inverse laplace transform of Y(s) 


32 
s(s +4)(s + 8) Poles: 


pl=0 p2=-4 p3=-8 


Y (s) = 


_ Ky K2 K3 
E DE) S 16447648) 


32 


"OTS, "T 


$——4 


K2 = -2 


l 3 
S sís + 8) 


32 
Ka = s(s +4) E. = 


1 2 1 
aant G44) 038) 


y(t) = (1 — 2e7*' + e79')u(r) 


PT ts ee TVA A T: 
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MATLAB 

cle 

s™s s 

num=3 2| 

den= poly([o di -8]) Yu as it ' Portal Fine 
[r.b,k]=residue (num, den) —— A cents) 

g= r(1)/(s-p(1)) + r (2)/(8-p(2))+ r (3)/(s-p(3)) 

pretty (91 

L Lap lace (a 


num = 
de 
den z 
1 12 32 0 
c = 
1 
-2 
1 
p a 
-8 
-4 
R 
k = 
LJ 
1 2 


ans = 
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Case 2. Partial- fraction expansion when F(s 


rey a NOS) 
ast = D(s) 
l Nis) 
~ (s+ py) G- p3): (5t Pa) 
Ky Kz K: 
= ————— qUe x ur ns 
(p (s any (s + p1) 
Kra dias K, 
(s + pa) (5+ Pn) 








Fi(s) = (s + p4) F(s) 


) involves Real and Repe 





ated poles 


The general expression for K1 through Kr for the multiple roots is 


1 d’'F,(s) 


g "cT 


ara "ipe 


ppt 


Ex: Find the inverse laplace transform of F(s) 


2 


FG) 7 e+e 


2 Ki K^ K3 


a G+) Ga2) T 642) 


K1 = (s + 1).F(s) = Eu em d 





K2 = (s + 2)2. F(s) = ca = —2 


K3 


4 407. a2 eA. 
L9 Fm =2 





f(t) =2€* — 2t! — 257? 


0! 1 














P d 


syms s 
F= 2/((s41)*(s42)^2); 
disp('F(s)') 
pretty (F} 

dispri'£itj') 
ilaplace(F) 


(S * 1) (s * 2) 


f(t) = 2/exp(t) - 2/exp(2*t) - (2*t) /exp(2*t) 
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Ex: Find the 1 inverse laplace transform of F(s) 
c^ 4 2s F4 


FS) > ~——= 
(s + py 


His b 
Fis) = 32. —1 e de 


de Der dis pns e 
MS) stl (s-ly (s+ 1) 


A(s) 
i t 2S » 3). d 


d + B(S) 
le (s + 1) ml i 


e 


dI + i 


d 


O 
T 


L 


2s + sx) 


ET 


Es * 1 M 


A 
ji 





Ro = L'S) 


| | oF U A | 
== y l — 4+ «P ! 7 * d 9 
t : + ] lc T s la T ni 


zo. Pe! 








—(I-15)e' ort <U 
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: i imaginary poles 
Case 3: Partial- fraction expansion when E(s) involves complex or imaginary | 


Ex: Find the inverse laplace transform of F(s) 


s$ + 2s 2 — (s +14 fl)(s + 1 — jl) 


1 Qi ais + ay 


Els) = —————-—4—————— 
(s) sL +2s +2) s .s?+25+2 


A = P EEES] +(ArS ras yS 





= ( S94 LL G ED 


un 


SC.s mas 42) 








EAS OTON) 


ME E | 
, Airaa sS A (201-435: S AZAL 


S C sn LS +2) 


ER 0» = € à 2 E 








“5 œ ¡e 
Oy » ab 
Ave des x = "mat 
l L 
A VERG Lud 
\ eR ls: Sd Motu Me rar OR 
> (S425 +2.) S Sa 258 AZ S * S ue ET 
S DA S Ar O 


2 e (1-4 


st | at 





SS 


A.J.A 
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Ex: Find the inverse laplace transform of F(s) 
$ + 2s? + 6547 











rs a TTS a spate] TTE 
Feces) L A is m "E A 
S a SAS” £ 14S45+2 œ L 
Fa ( s) I s EN MM - —À OL [s L 
r 2 c : 
(Sr KT + L (544 540% i 
. a [va [2 
s Fils). 43 4 me 2 «(t5 l2) 
( Sa X ) I l2. ) 
Por SL de HA. o Maleae 
i > 4 2 [Va "de 
[19 CS at. ) CELS 2 
a B dass A 
- - PG» : Y | ı CS) | RR E =5 
SE 
, iie 
DU ) es up, 48. us | 
= — st C ol 
dè Ma SYTA y 
I A DO) « Ste t 25% | 
la D Ua 
clc 
syms 3 


gs» s + 1 *( 2/(s^243845)); 

gt= ilaplace (gs): 

gt=simple (gt): 

C 

’ / 1/2 \ 


1/2 / E Y | 19 C 
4 19 exp] - | sin| ------- | 
2 / \ 2 i 
q(t)= dirac(t) + dirac(t, 1) + wr wren wn nn owe e e e e r m n a emm 
19 


x MATLAS fan chos l eat (S idle ; d oc | 
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—————————————— — — 


n——————————————— 





cA 
" Pu MZ — Y 
fa » Fr P E hA vesc. jrnrlac oe A bori 




















EsSu S A SAAL — AS -t = 
; SC s 1) 
= SL ua DA e o CMM. LA 
wl 
= E ZS Ss im ES i ETR 
$.— o m— S PL . 
Az - S4S E. sell EE Ayuso 
SS c korr. PI E TL E 
$--I So AS Tr Lio 
2$ 4.2 m idi. 
a Ce 
l 2 
d di. o 2 SA EX e 
| S S4 | 
bed. - SCO x al SU) + 28604 $ 7 3 e E eZ 
Lo: 
L K, Pf. el rasiolací mio, d Le) 
e) e i 
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Laplace Transform Solution of a Differential Equation 


EX: Given the following different; 
18 CMierential equation, solve for v(t) if all init; b 
Use the Laplace transform. d ; y(t) if all initial conditions are zero. 


dy dy 
TE + 12d + 32y = 32u(t) 


S?Y(s) + SY(0)- Y(0) + 12 ( SY(s) - Y(0) + 32Y(s) =32/s 


if all initial conditions are zero 


(S^- 128 + 32)Y(s) =32/s 





32 32 
A A A 
s(s* + -12s +32) sí(s+4)(s + 8) 
3 K, Kj», K 
TG) = s(s--Ay(s +8) $ Tn * (548) 
K = 1 


7-468). 


m 
2 ss +8), 


— pants 
e. 


3 
7 s(s +4) 


Lond 
— 


fa -B 





K3 


1 l 
¥(9) =; G44) * 48) 


y(t) = (1-20  e"*)u(r) 


MATLAB 


syms Ss TR Y(s) E^ M ins ----------- 
Y= 32/(s*(s*2 -t12*s432)): | J 
disp('Y(s)=') 3 (S + 12 s + 32) 
pretty (Y) 
yt=ilaplace (Y); 
disp('y(tj") 


y(t) = 1+ exp(-8 t) - 2 exp(-4 t) 
pretty(Gt) 


- , ` 
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.qnsform 
«y the Laplace trans 
le Solve using the 
Examp | T P : 
g YA 7e 


Application of the Laplace transform leads to l 


sY (s) — y(0) — Y(s) = 523 ; 
therefore. v ‘ i 1 1 ) 
2 clem a uem US —  —ÀÀ ——Ó . 
9\s5—1 s—3 


Y(s = 371 * (s— D(8-3) s—1 4 
Using the table to find the inverse Laplace transform, We obtain 
3 1 . 


y(t) ^ L a py] 226 (e. 5 
MATLAB 
cic 
syms y 
M us REP ae"? 


yt gm]: Dy-y= 
disp('y(t)=") 
pretty (yt) 

ezplot (yt, [- -10 10]) 
grid 

xlabel('t') 
ylabel('y') 


exp(3 ty2 + (3 exp(t)y2 
LIRE iso a 


U 
cot aaaea da esec e duos- 
a da g a =a RL coos gwh e p p. eo 
sb aaee de... < 
---- - da 


^ 
b owe oli Seeta AA A sna esie daa 
cn... -.¿4- Y ess 








A.J.A 
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Example . Solve 


mo, í : | 
y —3y -2y—6€*. y(0)=1, y(0) — 0. 
Applying the Laplace transform, I find 


2 
L" = 38 + 2)Y — s 3=- l í s? — 68 +10 
+32 3 > YO =D 


A b B " C 
ATA dl 











Y(s) - - 


To find A, B and C here is especially simple. tara example, for A multiply both sides by s— 3 and 
plug s = 3 into the expressions to obtain A — b. In a similar way B = —2 and C = 3 
Therefore, using the linearity of the inverse n trausform, we will find 


y(t) =v} =2 ze e, 


Example . Solve 
y" — 10y' +9y — BL y(0) = —1, y'(0) = 2. 


Applying the Laplace transform to both side, we find 


E 4.7 2 3 
2 ee A 5 + 12s* — 3 
(s^ — 10s + 9)Y +5-2-10= _ kai Y (s3) = 3iis - BY — 1) BI 


To find the inverse Laplace transform we will need first simplify the expression for Y (s) using the 
partial fraction decomposition: 











5 + 12s? — s? A B C D 
Sa x: Pee eo) 
We find 
5 31 50 
B=-,D=-2,C=—,A=-s- 
| 9’ imi 81 
Therefore, using the linearity of the inverse Laplace transform, 
50 5t 31 y, e 
a is Xu mU 
y(t) 81 +3 ts? e 
n pr. PES EE, 


syms Y r 
yt= -dsolve(' D2y-10*Dy+9*y=5*t",'Dy(0)=2"," y(0)=-1'); 
disp('y(t)=") 


OS L, A A EH 
«miissen 4 air dada dad: tabiat 3. iiaia AS oe vh "dcs 





pretty LYLE) 
ezplot (yt, [-10 10]) eee eHeemMeHbeH eR TD T 
grid ME EE e. po "» A me wmm cu qa... 
xlabel('t') w a d, y an A 5 L " ' 
nim d Wis se MN g at m a: 

x E GRE X D $13 4 €. 
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Example . Solve ML 
y” — Gy + 15y = 201038, y(0) = —1. y (0) = l. 


We have 





, | G 42099424 — As B ot. 
(só — 6s + 15)Y TR ET: => TI = (ay ETT s?--9 s^ — 6s +10 


To find the constants. we need to simplify the expression on the right (find the common denominator) 
and equate the coefficients at the equal powers: 

$:A-CLI-l 

s: -GA+B+D=2 

3:154 — 6B 4- 9C = —9 

s9: 15B +9D= 24 


The solution is 





1 11 5 
A= 19°18 =p =p =F 
Hence, we got 
xa A fayi. discos 
Visi = l mer T, 
(3) = 39 (5 Y Aus) 


Now we need to find the inverse Laplace transform. Let us start with the first term: 


s +1 S 1 S 1 3 ] 
ef stab aot Sm bibe s coasts! aa 


The second term is slightly more involved. Rearrange the expression in the following way (remember 


that we can always add and subtract the same expression and multiply and divide by the same 
expression different. from zero): 





—lis-25 — —1154+25 
|| —11(s—3)-3 
~ (s — 3)? +6 


(s — 3) 8 V6 


= 1 NO 
(s—3)?+6 V6(s—3)?4+6 


Now 


d —lis + 25 . 8 
6 I —————— = —11 3t v t — A 3t S] t. 
L E e cos V/6 Te sin V6 


The final answer hence is 


! 1 1 V 
y(t) gy) e 15 (cos 3t + 3 sin 3t — 11e% cos V6t — JA sin vit) | 
2 D 
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Week 72 (cont) 


V Mine 3 | . 
iodelinz in the Frequency Domain ......... ............EE362 


Mathematical model: Re 


system presentation of the input- output (signal) relation of a physical 


— — 
fu — 
EES 
















Sensor 









4. Implemenation 


Controller 


SS = AS D 
3. Design 


Physical Output 
in i 


" Modeling 





Mathematical model 


2. Analysis 


Remarks 
- A model is used for the analysis and design of control systems. 
- Modeling is the most important and difficult task in control system design. 
- No mathematical model exactly represents a physical system. 
- Do no: confuse models with physica! systems. 


- In this course, we may use the term system to mean a mathematical model. 


Math model Æ Physical system 
Math model x Physical system 
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Transfer function: A transfer function is defined by 


4 
G(s) i= r x 
U ( s) G(s) Y (s) 


The transfer function of a Linear, Time Invariant LTI, differential equation system 1s — 
as the ration of the Laplace transform of the output to the Laplace transform of the input under 
the assumption that all initial condition are zero. 








Let us begin by writing a general nth-order, linear, time-invariant differential equation, 





d'c(t) . dies) d" r(t) dare) 
Onn IO ll ++ + agc(t) = Do dm + Burl + «+++ bor(t) 


Where c(t) is the output, r(t) is the input, and the a;'s, b;'s, and the form of the differential - 
equation represent the system. Taking the Laplace transform of both sides, 


an C(s) + an-s”! C(s) +--+ + a0C(s) + initial condition 
terms involving c(t) 
= bms” R(s) + ba 457^ R(s) +--+ + bgR(s) + initial condition 
terms involving r(t) 
If we assume that all initial conditions are zero, 


(an5" + an 1577! + -+ + ag)C(s) = (bus + by as" a + +. + ba) R(s) 


Now form the ratio of the output transform, C(s), divided by the input transform, R(s): 





_ ins” + by a7! +--+ + bg) 
(a,5" + a, 507! +- + ao) 


ns anaa A  —— a a —— ee 





We call this ratio, G(s), the Transfer Function and evaluate it with zero initial conditions. 
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Definition: 
Poles of G(s) = roots of the denominator of G(s) 


= roots of characteristic equation 


Leros of G(s)= roots of the numerator of G(s) 


Stable: A system G(s)is said to be stable if the output c(t) of the system with a bounded 


input r(t) is bounded. 


The transfer function can be represented as a block diagram, as shown in Figure below, with 
the input on the left, the output on the right, the system transfer function inside the block. 






R(s) (b,,5 + bys! +--+ + by) 
(as a, 57 + +++ + ag) 


A block diagram of a system is a pictorial representation of the functions performed by each 
component and of the follow signals. 


Summary 


To derive the transfer function, we proceed according to the following steps: 
1- Write the differential equation for the system. 
2- Take the Laplace transform of the differential equation, assuming all initial condition are zero. 


3- Take the ratio of the output C(s) to the input R(s). This ratio is the transfer function. 


We can find the output, C(s) by using 
C(s) = G(s) . R(s) 


C(t) = inverse lablace (C(s)) = Response of the system to the input r(t). 


If input = step function output (step response) = c(t), 
If input = delta function output (impulse response) = c(t) , 
If input = ramp function output (ramp response) = c(t) 
A.J.A Fall 2017 
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Some standard forms for transfer function 


I- Ratio of polynomials 
m 


Md = bet St ot má 
Ga t C, S4 ha up Ars” 


2- Factored forms or product form 
y 29005 
| ($-39(5-29-.- (S-2n) 


K E= PAISES piin CS- Pn) 


MIS) < kK 


3- Come form : comes from the partial fraction expansion 


] | n - 
QI < kt == 
der S- 5h, 


MATLAB: Use MATLAB to generate the transfer function: 





_ S(s + 15)(s + 26)(s + 72) 


GLS) - . hei i oll he MN 
ts) sí(s + 55) (5? +27s + 52) 
In the following wavs: 
a- The ratio of factors b- The ratio of polynomials 


'Factored' 
Gp- zpk([ -15 -26 -72],[0 -55 roots([ 1 27 52])'],5) 
'polynomials' 


G= tf (Gp) 
Facored 
polynomials 
Zero/pole/gain:| 
5 (s+15) (s*26) (s+72) Transfer function: 


A e 5 S^3 + 565 s^2 + 16710 s + 140400 
s (st55) (st24.91) (s42.087) | | | 1 nn eee MMMM 


s^4 + 02 ad + 1537 s*2 + 2860 s 


———— — ————À —M — —À MM a P 
g ——À À— 
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Ex Use unction ; 
MATLAB to generate the transfer function in the following form (factors form) 


3 2 
BC) a es S" + 9s* -23s +15 


st + 2083 + 1405? + 400s + 384 





ratio form Gris) = 


clc 
num- [1 9 23 15]; Transfer function: 
den-[ 1 20 140 400 384]; s^3 +9 s^2 + 23 s + 15 


disp(' ratio form Gr(s) =') MEG uH AS + 140 s^2 + 400 s + 384 
Gr- tf (num,den) $ ratio 


[z,p,k]z tf2zp (num,den) Xa p k = 
disp(' factors form Gf (s) =') T 
Gf- zpk(z, -5.0000 -8. 1 
Pk(z,p,k) % facot 48 Gant -6.0000 
-1.0000 -4.0000 
-2.0000 
factors form Gf(s) = 
Zero/pole/gain: 
(3+5) (s43) (3+1) 
(5+8) (3+6) (s+4) (=+2) 
Ex Use MATLAB to generate the transfer function in the following form ( ratio form) 
G(s) (s+ 1)(s + 3)(s + 5) 
Sl = ————————————— 
(s + Z)(s + 4)(s + 6)(s + 8) 
num- poly([-1 -3 -5]); 
den- poly([-2 -4 -6 -8]) 
disp(' ratio form Gr(s) =') 
Gr= tf(num,den) % ratio 
Transfer function: 
S^3 +9 s^2 + 23 Ss + 15 
S^4 + 20 s^3 + 140 s^2 + 400 s + 384 
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MATLAB 5 
on BCs) /A(s) using MAT- 

Given the zero(s), pole(s), and gain K of B(s)/A(s), obtain the function B (s)/A(s) E 

LAB. Consider the three cases below. 

(1) There is no zero. Poles are at —1 + 2j and —1 ~ 2j. K = 10. 

(2) AzeroisatO. Polesareat -1 + 2/and -1 — 2j. ^ = 10. 


(3) A zeroisat -1. Polesareat -2,-4and -8. K = l7. 


z= [] e Pole-Zero Map 
, ECHO. MM 

p-[-1*2j ; -1-23]; | 

k=10; 

[num,den]- Transfer function: 

zp2tf(z,p,k); 10 2 

Gs tE (pum den) 3333  -———— i 

pzmap (num, den) wia dE se = 

z=[0]; 

z[-1423j3.; =1+23] 7 

ae 3 ji Transfer function: 

[num, den] = 10 s : 

SOIE A 7... .".: =. 4... E 

G= tf (num, den) 32 dE 8 d E 


pzmap (num, den) 








z=[-1]; 

p=[-2 ; -4 ; -8]; : B NM 
k=12; | | 

[num,den]- Transfer function: T 

zp2tf (z,p,k) ; 12 s + 12 s 

G=tE Lam - den) - == mesa K ox - y 
pzmap (num, den) 53 + 14 s^2 + 56 s +64 È 








———— ——— 9— a ae 








8 -6 4 2 
: Real Axis 
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Ex 


Find the transfer function represented by 


^f 
ací 


Y t 2c(t) — r(t) 


Taking the Laplace transform of both sides, assuming zero initial conditions, we have 


sC(s) + 2C(s) = R(s) 
The transfer function, G(s), is 
CS 1 
um) = (s) s+2 
Ex 


Use the result in EX1 to find the response, c(t) to inputs, assuming zero initial conditions. 
a- r(t) = u(t), unit step b- r(t) = delta(t), impulse c- LU - t, ramp 


a- LU 7 step function R(s) = l/s 
l 
C(s) = R(s)G(s) = sr 
1 1 1 


de s+2s  s(s-2) 


Expanding by partial fractions, we get 
1/2 1/2 


s+2 
Finally, taking the inverse Laplace transform of each term yields 





ett) = ; - =e 
b- r(t)- impulse function R(s) = 1 


C(s) = G(s).R(s) 





1 1 
C(s) = =p RG) aic 
c(t) = exp(-2t) 
c- r(t)=t,ramp function R(s) = 1/s? 
i d 1 
C(s) = au i a 
(5) s+2 s* s*(s + 2) 
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Expanding by partial fractions, we get 


1/4 T 1/4 


EIS s ..85 


c(t) »- 0.25 -- 0.5t +0.25exp(-2t) 


G= tf(1,[1 2]) 











step (G) 
figure Step Response 
- 0.5 T rM tM —ÀMRÓÓMMM ——— enne EE | 
impulse (G) AA | 
04 | 
| 
0.35 
0.3 
3 0.25 
E 
« 
0.2 
015] / 
| 
ie | 
0.05 | 
La "-———— e 
0 0.5 1 1.5 2 2.5 3 
Time (sec) 
Impulse Response 
1 NM EU ii TIN ÉÉ m q SOS era e = | 
0.9 
0.8 
0.7 
0.6 
2 0.5 
E 
T 
0.4 ; 
0.3 | 
0.2 | 
| 
K 4 
0.1 l T e j 
o 
Qi ————— e e NI 
0 0.5 1 15 2 2.5 
Time (sec) 
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Ex 
Find the transfer funct; ES 
unction, 9G) = C(s)/R(s), corresponding to the differential equation 
d*c La d € d dc t dtr dr 
d? de h a TA 


Taking the Laplace transform of the differential 
yields: 


equation assuming zero initial conditions 


in zi . P x 
s C(s) + 3s^C(s) + 7sC(s) + 5C(s) = S'R(s) + 4sR(s) + 3R(s) 


Collecting terms, 

— (s? +35 -- 7s +5)C(s) = (s^ + 4s + 3) R(S) 
C(s) — s'*4s*3 — 
R(s) s +3s°+7s5+5 

Ex 


Use the result in EX3 to find the response, c(t) to an input, r(t) = delta(t), a impulse function, 
assuming zero initial conditions. 


To solve the problem, we use the equation 


C(s) = G(s).R(s) 


Where 
C S +45 +3 
C + G 
R(s) s +35 +7s+5 
C(s) =R(s) . G(s) 
Since r(t) - u(t) step function , R(s) = 1, 
CS) s? --4s 4-3 " 
$) = —— > 
5343524 7s 4 5 (s) 
_ s24+45+3 (s) = Bis pecs 
(s34+3s2+7s+5) 53435247545" 


c(t) = inverse laplace of C(s) 
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i > C(s). 
Note: Using a hand calculation, ‘ind the inverse Laplace transform of C(s) 


Program 


cla 
syms s 
G= (s*2+4*54+3)/(8"3 +3%*s*2 +7*s+5) ; 3 2 

disp('G(s) = ') 

pretty (G) lake T PE 


R= 1 ; % where r(t) = delta function 

C= R * G; 

disp('C(s)= ') 

pretty (C) - 
disp('c(t)= ') R S 

c= ilaplace(C) ; S +33 +7 8s L B 
pretty (c) 


G [ET m 


expí(-t) (cos(2 t) + Sin(2 t)) 


Program 


G- tf(I1 4 3],[1 3 7 5]) 
impulse (G) 
Impulse Response 


me ——— o « — NN 
—— o ; — e e — 





Amplitude 





-0.4 lc qos om t mte € 
0 1 2 3 4 5 


Time (sec) 
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Ex5 
Repeat the EX3 if the input r(t) = u(t), step function 


C(s) = G(s).R(s) 


Where 
eu aaa a 
R(s) so +3s°+7s+5 — G(s) 


C(s) =R(s) . G(s) 


Since r(t)=u(t) step function , R(s)= l/s, 





2 
si+4s+3 
C(s) = ————————.R(s 
(s) s3+ 3544+ 75s4+5 (s) 
s(s3+3524+75+5) 53435247545" S 


c(t) = inverse laplace of C(s) 


Note: Using a hand calculation, find the inverse Laplace transform of C(s). 


Program 


S +4s + 3 
-F litt ie RL MEME 
syms s 
G= (s*2+4*st+3)/(s*3 *3*s^2 +7*s+5) ; see a 2 
disp('G(s) = ') 
pretty (G) 


R= 1/s ; % where r(t)= step function 2 

C= R * G; 

displ'C(s)9 ') 0 
pretty (C) 3 2 


disp('c(t)= ') 
c= ilaplace(C) ; c(t)= 


pretty LG) 
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Example: MATLAB Step Response | ETE 
Use MATLAB to find the step response of the following trans 
79s? + 916s 4- 1000 
Pal = Op 


; ials as such: 
We can separate out our numerator and denominator polynomial 


[79 916 1000]; 


num 


den = [1 30 300 1000 0]; $den- poly([0 -10 -10 -10]) 
BYE = Ef (num; den) 
T = 0:0.001:10 Step Response 


12 O TO O —— ^ 
‘ N ^ 


Step(sys, T) 

hold on 

plot (T, heaviside(T),'b') 
gtext ('input') 
gtext (’ output’ ) 





Time (sec) 


num = [79 916 1000); 
den = [1 30 300 1000 0]; den= poly([0 -10 -10 -10]) 
£ = 050.01:10; Simulation time = 10 seconds 


U = ones (size(T)); 

SYS = tf (num, den); 

lsim(sys, U, T) 

output) 

title('Response to unit step') 


input unit step 
construct a system model 
Simulate and plot the response (the 


o? o? oo o9 oo 


Response to Unit Step 
U 





' U 
I A p MEE E a o E E PT x 





Ame Hude 


E T ona a C crea MS T H NE 
LU 


Time (sec) 
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T 0:0.01:10» 


U = Chee size ps 
num = 









9. 
0 S i 
m «1/3 Big 1000] y jy ation time - 
are 7 11 30 300 1000 “Aput unit Ste ? Second; 
LOT Clio des d 3 
Y= lsim(sys, y a cs POLY LTO iy 70 
PLOE (E + D. Mages r | o Simulate “107 ) 
Citle('p Pr Y, ht and plot the r 
L ““PONSe to Unit "TRS (the output) 
“Send ("unit Step' : Step!) 
qrid r Output’) 
| a T |: e] 
: | ~~ Output 
wa 
7 8 9 10 
Ex 


For each of the following transfer function, write the corresponding differential equation. 


X(s) | s+2 
F(s) 33+ 11s? + 125418 


Cross multiplying, (s3-+8s2+9s+1 3)X(s) = (s+2)F(s) 


Px o dx de dW 
Taking the inverse Laplace transform, —7 + M + y^ + 15x= up -2Rt) 
| H | 
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